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ABSTRACT 

We numerically investigate light curves (LCs) of shocked circumstellar shells which 
are suggested to reproduce the observed LC of superluminous SN 2006gy analytically. 
In the previous analytical model, the effects of the recombination and the bolometric 
correction on LCs are not taken into account. To see the effects, we perform numerical 
radiation hydrodynamic calculations of shocked shells by using STELLA, which can 
numerically treat multigroup radiation transfer with realistic opacities. We show that 
the effects of the recombination and the bolometric correction are significant and 
the analytical model should be compare to the bolometric LC instead of a single 
band LC. We find that shocked circumstellar shells have a rapid LC decline initially 
because of the adiabatic expansion rather than the luminosity increase and the shocked 
shells fail to explain the LC properties of SN 2006gy. However, our synthetic LCs are 
qualitatively similar to those of superluminous SN 2003ma and SN 1988Z and they 
may be related to shocked circumstellar shells. 
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1 INTRODUCTION 



Origins of superluminous superno vae (SLSNe) wh ose exis- 
tence are recognized recently (see iGal-Yaml |2012| for a re- 
view) are still mystery. Several mechanisms to explain their 
huge luminosities are suggested so far fe.g.. ISmith et af 



20071 ; IWooslev. Blinnikov fc Hegerll2007l; iKasen fc Bildsten 



20101 ; Moriva et all 120121; ICha tzopoulos, Wheele r, fc Vinko 
20121 ; iDexter fc Kasenl |2012| ; lOuved et ail I2012T ). At least 
for a subclass of SLSNe which is named as 'SLSN-II' in 
IGal-Yaml (|2012T ). ejecta-circumstellar medium (CSM) inter- 
action is likely to be a main power source because most of 
them show narrow emission lines which are supposed to be 
rela ted to the existen ce of the dense circumstellar shell (see, 
e.g.. lFilippenko1ll997h . SNe can be very bright at their early 
epochs with the dense circumstellar shell since the kinetic 
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energy of SN ejecta can be efficiently converted to radiation 
due to the interaction. 

Dense circumstellar shell causing SLSNe are likely to 
be d ense enough to cau se the shock breakout in the shell 
(e.g., iMoriva et ai1l2012l ). There are two possible locations 
at which shock breakout can occ ur in a dense shell. The tw o 
cases are clearly summarized in IChevalier fc Irwinl l|201ll) . 
In one case, the optical depth of the dense shell is high 
enough to keep photons in the forward shock until it gets 
near the surface of the dense shell. In other words, the shock 
breakout occurs near the surface of the dense shell. On the 
other hand, if the optical depth is low, the shock breakout 
can occur within the shell well before the shock reaches the 
vicinity of the shell surface. In this case, photons emitted 
from the shock diffuse in the remaining shell after the shock 
breakout. LCs from the latter cases have been numerically 
investigated and t hey are found to be consistent with the 
LC o f SN 2006gy l|Moriva et al.ll2012l ; ICinzburg fc Balberd 
120121 ). 
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Table 1. List of initial conditions 



NellTIC 


km s _1 


Mq 


10 15 cm 


rp d 

10 4 K 


composition 


M01 


4,000 


10 


2.4 


1 


solar 


M02 


4,000 


10 


2.4 


4 


solar 


M03 


4,000 


20 


7.2 


1.7 


solar 


M04 


2,000 


20 


7.2 


1.7 


solar 


M05 


8,000 


20 


7.2 


1.7 


solar 


M06 


4,000 


20 


0.72 


13 


solar 


M07 


4,000 


20 


7.2 


1.7 


C 50 % + O 50 % 



a velocity of the outermost layer of the shell 

" mass of the shell 

c radius of the shell 

^ initial temperature of the shell 



ISmith fc McCravl (|2007h (SM07 hereafter) investigate a 
SLSN LC model resulting from a circumstellar shell through 
which a shock wave has gone through. This model basically 
corresponds to the former picture of the shock breakout, 
i.e., shock breakout at the surface of a dense shell. In other 
words, SM07 consider a dense shell from which photons start 
to emit after the passage of a shock wave. They apply the 
LC model of adiabatically cooling SN ejecta formulated by 
lArnettl < |l98Cf >. By simply comparing the shape of the model 
LC to the 7?-band LC of SN 2006gy, SM07 concluded that 
SN 2006gy can result from a shocked circumstellar shell. 

However, there are many simplifications involved in 
the SM07 model which are not discussed so far. The 
existence of the recombination wave in shocked H-rich 
CSM is presumed to affect the LC as in the case of 
Type I IP SNe (e.g.. ICrassberg. Imshennik. fc Nadyozhinl 



19711; iFalk fc ArnettJ Il977l; iKasen fc Wooslevl 120091 ; 
Bersten. Benvenuto. fc Hamuvl 1201 ll ). In addition, SM07 



compare a bolometric LC obtained from their model to the 
fl-band LC of SN 2006gy. As the shocked shell should have 
temperature close to 10 4 K at t he beginning to ex plain the 
observed SN 2006gy properties l|Smith et al.ll201oh . a large 
fraction of emitted photons is not in the R band and the 
bolometric correction should be considered. 

To see the importance of these neglected effects in the 
SM07 analytic model, we numerically follow the system 
suggested to explain SN 2006gy in SM07 by using a nu- 
merical radiation hydrodynamics code STELLA. STELLA is a 
one-dimensional radiation hydrodynamics code. STELLA can 
treat multi-frequency radiation and the opacity is calcu- 
lated based on the physical parameters under the assump- 
tion of the local thermodynamic equilibrium. This assump- 
tion is a good approximation to obtain a LC from a shocked 
shell because of its large density and optical depth. Thus, 
STELLA is a suitable code to see the effect of the bolometric 
correction and recombination which are not taken into ac- 
count in the SM07 analytic mo del. For the details of STELLA , 
we refer to other articles (e.g., Blinnikov fc Bartunov 1 19931 ; 



iBlinnikov et al. 2000, 2006; Blinn ikov fc Tolstovll201ll 



We start by showing the initial conditions of the models 
which are constructed based on SM07 in Section [2] and then 
show the numerical results in Section [3] Discussion is pre- 
sented in Section [4] and we conclude in Section [5] We apply 



the same distance to the host galaxy (73.1 Mpc) and extinc- 
ti ons (Galactic Ar_= 0.43 mag + host Ar — 1.25 mag) as 
in lSmith et all (|2007l ) to the observed LC of SN 2006gy. 



2 MODELS 

Initial conditions of our numerical calculations are con- 
structed based on SM07 at first. We do not follow shock 
propagation in shells to make the initial conditions the same 
as those of SM07. The initial conditions are supposed to re- 
sult from the shock passage in the dense circumstellar shell. 
Since we do not treat the shock wave, the smearing term 
in STELLA code which affects the conve rsion efficiency from 
kinetic energy to radiation discussed in lMoriva et al.l (|2012l ) 
does not affect the results obtained in this paper. 

Table[T]is the list of our initial conditions. The initial ra- 
dius of the SM07 model suggested for SN 2006gy is 2.4 x 10 15 
cm. The mass of the shocked shell in the model is 10 Mq. 
We assume that the system is homologously expanding and 
the outermost layer velocity is 4,000 km s _1 , as assumed 
in SM07. The initial temperature is set constant in the en- 
tire shell. We try two temperatures for the SM07 system, 
namely, 10 4 K (Model 01 or M01) and 4 x 10 4 K (M02). 
The composition is solar in these models. 

We also investigate several configurations other than 
those suggested in SM07. M03 has the same velocity but 
the initial radius is three times larger than that of the SM07 
model. The temperature is set to 1.7 x 10 4 K to match the 
observed luminosity of SN 2006gy. The mass is increased to 
20 Mq in M03 to keep the shell optically thick. M04 and M05 
have the same density structure as M03 but the velocity is 
0.5 and 2 times of M03, respectively. We also show results 
of M06, which is more compact than the SM07 model. The 
composition is solar in M03-M06. M07 has 50% of carbon 
and 50% of oxygen, while other properties are the same as 
in M03. 



3 SYNTHETIC LIGHT CURVES 

Figure [TJi shows the bolometric LCs of M01-M03. At first 
glance, we find that no LCs are consistent with that of 
SN 2006gy. Furthermore, R-b&nd LCs shown in Figure QJ) 
are found to be more different from the .R-band LC of SN 
2006gy. This indicates the importance of the bolometric cor- 
rection. We discuss the LC behaviors in this section but most 
of our discussion can be found in the previ o us studies, e.g., 
ICrassberg. Imshennik. fc Nadyozhinl i| 19711 ); IFalk fc Arnettl 
( 1977D- 

The bolometric LCs start with the initial peak. The 
peak bolometric luminosity is Air R^crT^x, where a is the 
Stephan-Boltzmann constant. At first, the bolometric lu- 
minosity decreases due to the adiabatic expansion of the 
shell. If we assume that the homologously-expanding shell 
is radiation-dominated at early phases, the bolometric lumi- 
nosity should decrease following oc t~ 2 . This rapid decrease 
in the bolometric luminosity appears in our numerical mod- 
els. 

SM07 suggest that the bolometric LC of shocked shells 
would rise following oc t 2 at first because shocked materi- 
als expand homologously and the shocked shell is just an 
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Figure 1. Bolomctric (a) and i?-band (b) LCs of shocked shells. 
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Figure 2. Evolution of the effective temperatures. 



Figure 3. Evolution of the photospheric radii. The arrows in the 
figure shows the directions of the evolution of the radius with a 
given constant velocity. 



expanding blackbody. However, if we take into account the 
decrease in the blackbody temperature of the shell due to 
the adiabatic expansion and the lack of any heat sources, 
the effect of the temperature decline on the bolometric lu- 
minosity (oc t~ 4 in radiation-dominated shells) is larger than 
the effect of the radius increase on the bolometric luminosity 
(oc t 2 ). In fact, our synthetic LCs do not show the luminosity 
increase and the luminosity just declines. 

After the initial rapid luminosity decline, bolometric 
LCs start to be affected by photons di ffused in the sh ell and 
begin to follow the diffusion model of lArnettJ (|l98fj ). From 
this point, the SM07 analytic model starts to work. M01, 
whose initial temperature (10 4 K) is close to the blackbody 
temperature of SN 2006gy at the LC peak, is too faint at 
this epoch to explain SN 2006gy because of the initial rapid 
luminosity decline due to the adiabatic expansion. With the 
initial configuration suggested by SM07, the temperature 
should be around 4 x 10 4 K (M02) to explain the luminosity 
of SN 2006gy but it is inconsistent with the observed black- 
body temperature (Figure [2]). M03 has a larger radius than 
those of M01 and M02. Thus, the required temperature to 



get the same luminosity is small (1.7 x 10 4 K) and it is close 
to the observed values. 

Although the bolometric LCs at these epochs seem to 
follow the observed bolometric LC, we should be careful be- 
cause the bolometric LC of SN 2006gy is obtained by the 
7?-band LC without the bolometric correction. We need to 
compare LCs in the R band (Figure[]J>) . We find that the nu- 
merical i?-band LCs do not match the observed R-h&nA LC 
even in the models which give a good fit in Figure This is 
simply because of high temperatures in the shell and most of 
the emitted photons are not in the i?-band. We note that the 
strong Ha line observed in SN 2006gy is in the i?-band and 
the direct comparison between our numerical J?-band LCs 
and the observed J?-band LC can be inappropriate. How- 
ever, the Ha luminosities of SN 2006 gy at the epochs w e are 
interested in is just ~ 10 erg s" - 1 jSmith et alJl201Ch and 
the bolometric correction remains to be important. 

Another unavoidable and important consequence of the 
SM07 model which is not discussed in SM07 is the existence 
of the recombination wave in the shocked shell. In the SM07 
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Figure 4. Bolometric (a) and i?-band (b) LCs of the models with different initial velocities. 



model, there are no energy sources in the shell because the 
shock has already passed the shell and the shocked shell just 
cools down. At one epoch, the temperature should reach 
the recombination temperature as is the case for Type IIP 
SNe. This is not the case for the continuous ejecta-CSM 
interaction models because there remains an energy source 
(shock waves) which can keep the shell ionized until the 
shock wave goes through the dense shell. 

The effect of the recombination can be seen by compar- 
ing M03 and M03co in Figure [T]i. M03 is calculated with 
our standard opacity table which takes recombination into 
account. M03co (M03 constant opacity) is calculated by forc- 
ing the scattering opacity of the system to be 0.34 cm 2 g~\ 
which corresponds to the fully ionized solar composition ma- 
terials. At first, when the shell is above the recombination 
temperature, the two LCs follow almost the same track. 
Then, two LCs start to deviate when the outermost layer 
reaches the recombination temperature at around 40 days 
since the LC peak (Figure |5J. The recombination wave, and 
thus the photosphere, move inside (in Lagrangian sense) af- 
ter this epoch. They eventually reach the bottom of the 
shell and the LC suddenly drops. On the other hand, the 
LC with the constant opacity continues to decline monoton- 
ically, roughly following the SM07 analytic model. Figure [3] 
shows the photospheric radii of the models and the effect of 
the recombination is clear. 

Another important consequence caused by the existence 
of the recombination is the strong dependency of LCs on the 
shell velocity. The epoch when the outermost layers reach 
the recombination temperature and the recession velocity of 
the recombination wave in the shell are affected by the shell 
velocity. This is simply because adiabatic cooling becomes 
more efficient in faster shells. M04 and M05 have slower and 
faster shells, respectively, than M03 and their LCs are pre- 
sented in Figure 2] At first, the LCs are expected to differ 
when the recombination start to play a role in the shells. 
However, the difference at this epoch is not significant ac- 
cording to our calculations. The time of the drop in the LCs 
clearly differs and the faster shells have earlier drops due to 
the faster recombination wave. 




days 



Figure 5. Bolometric and multicolor LCs from M06. The R- 
band LC has a rising phase similar to that of SN 2006gy because 
of the initial high temperature 1.3 X 10 5 K. However, the high 
temperature is inconsistent with the SN 2006gy observations. 



4 DISCUSSION 

4.1 Model with Initial R-band Luminosity 
Increase 

All the models we have presented so far do not have a phase 
with luminosity increase and the luminosity just declines. 
However, it is possible to have a rising phase in optical LCs 
from a shocked shell. Figure [S] shows one example of LC 
obtained from M06. The bolometric LC and optical LCs 
as well as the .R-band LC of SN 2006gy are shown in the 
same figure for the illustrative purpose. The evolution of 
the bolometric LC does not differ so much from the previous 
models but optical LCs of M06 have a rising phase. This is 
because of the initial small radius and high temperature. 
The optical luminosities are low at the beginning due to the 
initial high temperature. Then, as the adiabatic cooling is 
efficient due to the initial small shell radius, the shell cools 
quickly and optical luminosities increase accordingly. Then, 
the J?-band LC can be similar to that of SN 2006gy, although 
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days since the LC peak of SN 2006gy 
Figure 6. Bolometric (a) and i?-band (b) LCs from a C+O-ri 



the photospheric temperature is much higher in M06 and it 
is inconsistent with the SN 2006gy observations. 



4.2 Possible Corresponding Supernovae 

LCs of shocked shells obtained by our numerical calcula- 
tions have an initial rapid decline followed by a relatively 
long plateau. Although these features are not seen in SLSN 
2006gy, S LSN 2003ma (SL SN-II) qualitatively shows similar 
features (|Rest et al.ll201l[ ). The LC of SN 2003ma is differ- 
ent from other known SLSNe. The LC of SN 2003ma has a 
quick rise and quick decline followed by a long plateau phase 
which lasts for about 100 days while other SLSNe evolves 
more slowly. Then the LC drops by about 1 mag in opti- 
cal and the luminosity stays almost constant for about 1000 
days after the drop. The initial rise and decline as well as 
the plateau phase which lasts for about 100 days can be seen 
in some synthetic optical LCs obtained in this study (e.g., 
Figure [Tp), but the plateau phase after the drop which lasts 
for about 1000 days requires another emission mechanism 
like a continuous CSM interaction. 

SN 1988Z has a similar feature to SN 2003ma, al- 
though the luminos i ty is about 3 magnit ude smaller (e.g., 
iTuratto et"ai] Il993l ; lAretxaga et all ll999l L Because of the 
LC simil arity, SN 1988Z can also b e related to shocked shells 
(see also lChugai fe Danzi gcr 1994). Depending on, e.g., radii 
and temperatures of shocked shells, their luminosities can 
vary. There can be many other similar SNe with variety of 
luminosities and plateau durations, depending on the shell 
properties. 



4.3 Other Effects 

4-3.1 Shell of Carbon and Oxgen 

We examine a LC from a shocked shell with 50 % of car- 
bon and 50% of oxygen. A subclass of SLSNe is known to 
have no hydrogen features in their spectra and their com- 
positi on is likely to be d ominated by carbon and oxygen 
(e.g., iQuimbv et al1l201ll ). Figure [6] shows the results. As 
the opacity significantly decreases compared to the cases of 




days since the LC peak of SN 2006gy 
shocked shell (M07). M03 LCs are also shown for a reference. 



the solar composition, the LC declines much faster. In addi- 
tion, due to the high recombination temperature of carbon 
and oxygen, the effect of the recombination is less significant 
than the cases of the solar composition. However, we do see 
the significant drop in the LC due to recombination if we 
compare it to the LC model (M07) with a constant opac- 
ity (0.1 cm 2 g _1 ). Thus, we should also be cautious about 
the applicability of SM07 model when we try to apply it to 
carbon and oxygen dominated system. 

4.3.2 Opacity Table 

We also investigate the effect of the different opacity ta- 
bles adopted in STELLA (Figure[7]). The detailed descriptions 
about the three opacit y tables adopted here can be found in 
iTominaga et al.1 l |201ll ). 

In one opacity table ('inner shell' in Figure^, the inner- 
shell photoionization is additionally taken into account. The 
opacity table used so far in this paper assumes that all atoms 
and ions, except for hydrogen, are in the ground state to 
obtain the opacity from the bound-free absorptions as in 
lEastman fc Pintol |l993). In the other opacity table used in 
this section ('bound free' in Figure 0, exited levels are also 
taken into account in bound-free absorptions. 

As can be seen in Figure [7| in which LCs from the three 
opacity tables are shown, the difference caused by the dif- 
ferent opacity tables is small. Thus, the effect of inner-shell 
photoionizations and excited levels in bound-free absorp- 
tions on LCs is negligible in the system we are interested 
in this paper. 



5 CONCLUSIONS 

We numerically investigate LC properties of shocked shells 
which are suggested to account for SN 2006gy by SM07 an- 
alytically. We show that shocked shells fail to explain the 
rising part of the SN 2006gy LC because of the adiabatic 
cooling. In addition, in the declining part, we show the im- 
portance of the effects of the bolometric correction and the 
recombination which are not taken into account in SM07. 
SM07 compare their analytic bolometric LCs to the observed 
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Figure 7. Bolometric (a) and R-band (b) LCs with different opacity tabfes. 'M03' is computed with the default opacity table. The 
opacity table used in 'M03 inner shell' additionally takes into account the inner-shell photoionization. In 'M03 bound free', some exited 
levels are added in bound-free absorptions to the default opacity table. 



i?-band LC of SN 2006gy but the high temperature of the 
shell which is estimated from the spectral observations is 
against this assumption and i?-band LCs of shocked shells 
becomes flatter than bolometric LCs. Recombination also 
makes the LCs flatter than those analytically suggested in 
SM07. 

We also show the effect of the expansion speed of 
shocked shells. It affects the propagation velocity of the re- 
combination wave in the shells and changes duration of the 
LCs. The composition of the shells alters the opacity and is 
shown to change the LCs very much. 

Although we show that the LC of SN 2006gy is not con- 
sistent with the shocked circumstellar shell model, this kind 
of system can exist. We suggest that SLSN 2003ma and SN 
1 988Z may come from a shocked shell because of their qual- 
itative LC similarity to our numerical results. In addition, 
our numerical modeling indicates that shocked shells can be 
bright (sometimes more than ~ —23 mag in optical) and 
their luminosities can drop more than I magnitude in opti- 
cal within I day (M03). This kind of object can be observed 
in the future transient surveys and they can fill the bright 
and f ast declining pa rt in the explosive transient phase space 
(e.g., iKulkarn] 12012V The existence of such transients from 
a shocked shell indicates the existence of the explosive mass 
loss just before the explosions of some kinds of stars which 
is currently not known well and they can provide us a clue 
to understand such mass loss. 
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